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Whether the much studied X(3872) is an axial or tensor resonance makes an important difference 
to its interpretation. A recent paper by the BaBar collaboration [l[ raised the viable hypothesis 
that it might be a J^^ = 2 ^ state based on the Tr^Tr^Tr" spectrum in the X J/ip to decays. 
Furthermore, the Belle coUaboratioiipublished the tt+tt" invariant mass and spin-sensitive angular 
distributions in X ^ J/ip p decays Q]. 

Starting from a general parametrization of the decay amplitudes for the axial and tensor quantum 
numbers of the X, we re-analyze the whole set of available data. The level of agreement of the two 
spin hypotheses with data is interpreted with a rigorous statistical approach based on Monte Carlo 
simulations in order to be able to combine all the distributions regardless of their different levels of 
sensitivity to the spin of the X. 

Our analysis returns a probability of 5.5% and 0.1% for the agreement with data of the 1^+ 
and 2 hypotheses, respectively, once we combine the whole information (angular and mass 
distributions) from both channels. On the other hand, the separate analysis of J/V' p (angular and 
mass distributions) and J/^p cj (mass distribution) indicates that the 2 ^ assignment is excluded at 
the 99.9% C.L. by the former case, while the latter excludes at the same level the 1"^"'" hypothesis. 
There are therefore indications that the two decay modes behave in a different way. 

PACS: 13. 25. Ft, 14.40.Rt 



I. INTRODUCTION 

Although the X{3872) resonance is the most studied among the exotic XYZ states, since its discovery in 
2003, its quantum numbers have not been definitively identified yet. 

The CDF collaboration concluded from the analysis of the angular distributions and correlations of the X 
decay products that the possible quantum numbers are JP^ = 1++ or 2 '^[31. Similar results have been found 
very recently by the BELLE collaboration Q. In the latter paper the TT'^n~ invariant mass distribution is 
analyzed under the l"*""*" and 2 hypotheses, finding preference for the former whereas no preferred assignment 
emerges if an interfering contribution with the isospin-violating decay lj — >■ tt+tt" is added to the amplitude. 

The picture becomes more puzzling as one considers the analysis by the BaBar Collaboration of the decay 
X — > J/^jJ 7r+7r^7r° The expected Stt invariant mass distribution agrees with data slightly better if the 2 
signature is assumed. This result on the 3tt spectrum was later confirmed in 

In a recent paper Q the pion invariant masses in the decays X — ^ J/ip 2tt and X — )■ J/ip 3tt are simultaneously 
analyzed with a combined fit, concluding that present data favor the 1^^ assignment. In our view, the 27r 
invariant mass distribution is not able to resolve the two hypotheses despite of the high statistics, and the way 
the fit was performed leads to the dilution of the sensitivity of the Btt channel. 

To improve the analysis of all available data sensitive to the spin of the X i) we write the decay matrix elements 
for both the 1"''"'" and 2 hypotheses as given by enforcing Lorentz invariance and parity considerations, ii) we 
give a functional dependency on the decay momenta to the couplings introducing a length scale R - which is to be 
related to the finite size of the hadrons participating to the interactions, iii) we use the matrix elements of p and 
id decays to take into account the appropriate decay waves; we do not pursue the Blatt-Weisskopf description 
as we find that all spin structure can appropriately be taken into account with no further approximations, 
iv) we perform a global fit to exploit the information contained in all the distributions available and we adopt 
a statistical approach appropriate when distributions with different sensitivities to the parameters of interest 
are combined. 



II. MATRIX ELEMENTS 

The matrix elements describing the amplitudes X — > J/iJj V (where V = p,uj) are obtained by Lorentz, gauge 
invariance and parity considerations leading to the formulae reported below 
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In the X(l++) case we have 

{^Pie,p)Vi7^,q)\Xi\,P))^g,^v e'^-"- A^(P) eUp) r,;{q) (1) 

the polarization vectors carry a complex conjugation when referred to final states. 
In the 2 ^ case we have a more complicated structure 

{■ik{e,p)V{Ti,q)\X{^,P)) ^ g^^v Ta + g'2^v Tb (2) 

where tt is the polarization tensor for a spin two particle with mass"*^ and a standard notation is used for the 
remaining polarization vectors. We find that and Tg are given by 

TA^e'^ [p) (P) e^''^^ p.qpi^l{q)~{t,p^f^,q) (4) 

and 

Tb = Q" e'^'^P^ P, e;{p) rj^q) (5) 

where Q = p — q and P = p + q. 

The coupling giipv is real whereas the couplings (72i/>y and 521/)^ ^'"^ separately real but can have a complex 
relative phase. This is due to the fact that, on the basis of Lorentz invariance and parity conservation, we can 
indeed write three terms T^*^ , and T^*^ where A labels one out of the 3x3x5 polarization combinations 
which define Ta , Tb , Tc ■ Tq is the same as Ta but with a plus relative sign between the two terms on the rhs 
of (HJ. Only two out of these three terms are linearly independent (say Ta and Tb) for 



3x3x5 

Ey(>)y(>)- 
A C 



Ei^A^i'Ei^c^i' (6) 



A 



which is the equality of the Schwartz inequality: this holds if ziTa — Z2Tc where zi and Z2 are two complex 
numbers both different from zero. Thus we can exclude Tc and retain Tb and, in general, ziTa to characterize 
the decay amplitude. 

In Refs. [JHjIHI the P-wave fit functions contain a Blatt-Weisskopf angular momentum barrier factor of the 

]^ /2 

form (1 + R^q*^) , where q* is the X decay 3-momentum. The value of R cannot be extracted from data 
since the P-wave distribution will approach the S'-wave distribution in the limit P — > cx), so that if we let R 
free, the fit will not converge. On the other hand, in our discussion we do not need any barrier factor the decay 
wave being dictated by the expressions of the matrix elements. We instead take into account the finite size of 
the X (and of V and J/V' as well) introducing a 'polar' form factor, namely 

(7) 



(l-fP2g*2)n 

where g stands in general for gi^v, g2ipv and g'2^y- We tested the values n = 1 and n = 2 (the latter coinciding 
with the Fourier transform of a an exponential g{r) ^ exp(— r/P) strong charge distribution). Both the fitting 
functions turn out to be rather effective at improving our results, with no significative change for the the two 
choices of n. We also underscore that gi^v (regulating the S'-wave decay) is assumed to have the same polar 
behavior since Eq. ([7]) does not concern any orbital angular momentum considerations. The size parameters Rj 
will eventually be fitted from data. 

As for the p and uj decay amplitudes, we use 

{7r+ip)7:-{q)\pie,P))=gp2.e-p (8) 

which describes a P-wave decay (the square modulus of this matrix element is 5p27r times the decay momentum 
squared). For the lo we have 

{7r+{p)7T-{q)7r°{r)\u; (e,P)) = .9^3. e^"''''e^p,qpr, (9) 

The last two couplings will simply be written in terms of the partial widths r(p — > 27r) and T{uj — )■ 37r), as 
shown in the next section. 



^ The sum over polarizations is 

^7r^,^(fc)7r*^(A;) = -TpaT^/j + -T,,^T^q - -Tp^Tc/j (3) 

pol 

with T^i/ = —g^iu + k^ki^/m? and k'^ = m? . 
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III. DECAY WIDTHS 



We have to calculate the partial widths r{X J/ip tt+tt") and r{X J/ip tt+tt^tt^). In what follows we 
will neglect the p-ui mixing since we demonstrate in Appendix E] that it does not alter significantly the results. 
The partial widths in the narrow width approximation are 

^ pol , 

1 1 r 

.ydci>(2)^i(.+.-ip(.))p 



pol 



where by we mean the 2-body phase space measure. The decay momentum q* in the matrix element 

{tfj p{s)\X) coincides with q* = p*(to^, m^, s); similarly for the width in J/tp uj discussed below. 
The sum over polarizations in simply yields 

^ I (tt+tt- \pis)) p = p* (s, ml,mir (11) 

pol 

Finally, we can eliminate the coupling by evaluating the fjlip on the mass-shell and by relating it to the partial 
width r (p — ^ tttt) 

= 6ml Tip ^ nn) ^ (12) 

Inserting the above expressions (fTTj) . (|12p into (fTO|) gives 

r{X^J/i; ^+^') = 2jVi8W / dsJ2\{^Pp{s)\X)W{mjc,ml,s) 



pol 



(13) 



TT (s - m2)2 + {mpTp)^ y/s \p*{m'^p,ml,ml) ^ 
Similarly, for the uj we obtain 



^ " pol 

1 1 

X 



pol 



(14) 



Summing over the uj polarizations 



^ = [K + ^ - 27^0.) {uo' -ml- 4x2) _ (^2 _ ^2)] ^ ,23^^(7^) (15) 

pol 

where w = i?^o, x = ^ (Et,+ — i?7r-) and mo = TOttOj '^1+ = = m^-. An adimensional coupling has been 
formed by substituting (7^3^ — > -^g^^^ ^- We eliminate the coupling by evaluating on the mass-shell 

5^3^ ^ 6m^ r(c^ d$^'^M(m^)^ (16) 



This rescaling is arbitrary and in part relies on the narrowness of the ui. Avoiding the introduction of the term, the separate 
fit of Stt does not change significantly (for example, with n = 1, the 2 is unchanged, the 1++ gets worse from = 9-9 — ^ 11.1). 
In the combined fit of invariant mass distributions both hypotheses become a bit worse (again for n = 1, 2~ + : = 17.7 — > 18.4; 
x^ = 25.2 26.2) so that the Ax^ remains almost unchanged. See Sec. IIV Al 
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Inserting (fT5|) . (|16l) into (fT4)) we obtain 



^ pol 

1 m^^r^^ >57^(a; -> Stt) $^^^'(v^, mo, m+, 7n+) 



(17) 



TT (s - )2 + {mJi^^Y' Too, m+, to+) 

where 

$(^''(^5, mo, m+, m+) = / / [(m^ + s - 2^5 [u? - - 4a;^) - 4m^ [u? - m^)] 



mo 



with w„j = (ttIq — 4m^ + s)/(2y^) and 



(18) 



1 / (w2 _ (^^^^^ _ ^) , 



2Y 4m^ + (t^™ - tj) ^ ^ 

to be compared to the notations used in Q. 

In formulae and ((T7)l angular correlations are not taken into account for we factorize matrix elements. 
On the other hand, the only way to consider both off-shellness and angular correlations is to compute the full 
matrix element for the 1 — )■ 5 decay 



YiB^KX^KJI^P^ Kin- ~ / n (^(2.)^. (p. - E 



Vi 



(20) 



X 

pol 



The matrix element can be decomposed as 



(i^/+r 7r+7r-|S) = (/+r|^> — (7r+7r-|p> 



(21) 



X ^ ^i- — — 

Vx ~ "^x + *"^xl X 

We already gave a form to p\X) and (tt+tt^Ip) in Sec. [ll] for both 1++ and 2 ^ . Moreover 

(;+r|^) = (6^)^s(r)7^i;(;+) (22) 

(tt+tt^Ip) ep (23) 
el' ■ for 1" 

{i^*x)8^[VKf{pKr for 2- 



{XK\B) ^{ ^ '^'^ ^ ^ (24) 



The sum over polarizations of inner legs returns the usual numerators of spin-1 and spin-2 propagators. The 
expressions obtained are not reported here because of their algebraic complexity All matrix elements but the 
{ip p\X(2 '")) can be written in terms of one coupling only times a scalar function of momenta and polarizations. 

The widths are expressed in terms of sums of products of couplings times scalar functions of momenta. 
Products of couplings are absorbed within the fit parameters r j"^ see Sec. IIVI 

Since the integral in Eq. ([20]) has to be evaluated numerically, we use Monte Carlo techniques. The integration 
on the 5-body phase space is carried out using importance sampling on the Breit-Wigner peaks. As a further 
option of the calculation, unweighted decay configurations can be generated according to full matrix element 
weights. 

As usual, in the calculations we use the so-called comoving width, i.e. we rescale the width in the denominators 
of Eqs. (IT^ . p7|) and (|21l) according to the standard prescription mF — > {s/m)r. 



^ These are available upon request in form of Fortran routines. 
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IV. DATA RE- ANALYSIS 



In order to extract information on the spin of the X particle, we re-analyze the B — )■ XK data published in 
Ref. (for X J/V- TT+TT-) and in Ref. Q (for X ^ J/ip tt+tt-ttO). In particular, in the X J/iP tt+tt" 
sample we consider the di-pion invariant mass {m2-n) and the angles defined in the X rest frame in Ref. Q 
and described in Fig[T] the angle between the J/V' and the direction opposite to K (Ox), the angle between 
the 7r+ and the direction opposite to K (x), and the angle between the Z+ produced by the decay of the J/ip 
and the z-axis of a coordinate system where the x-axis is the direction opposite to the K and the y-axis is the 
component of 7r+ orthogonal to K {9i). In the X — )■ J/ if) 7r+7r~7r°, instead, only the three pions invariant mass 
is considered (ra^j,). 




With such distributions at hand the parameters of the model we have sensitivity upon are 

• the radii Rj (J — 1, 2) as defined in Eq. to which only the invariant mass distributions are sensitive 

• the relative amplitudes and phases of the two contributions in case J = 2 {g2ipv and g'2^y in Eq. ([2])). We 
redefine 

9 9 
92^,^ = cos ^, 32^5 = sin -j-e'f'^ (25) 

where ^ can be p, u or ang, depending on whether we are fitting the r7i27r, Tn^i-n, or angular distributions 
respectively. The parameters 9p and ipp satisfy 9p = ^ang and ipp — v^ang since in the p channel we have 
information on both mass spectrum and angular distributions. Only the angular distributions are sensitive 
to the 9 and if parameters 

• the overall normalizations, and ri^ (= giipp in Eq. ((HJ)). It has to be noted that the normalization 
depends on the distributions being studied 

Given the different sensitivities of the mass and angle distributions and the computational complexity of the 
angular fits, the combined fits are performed in three steps: i) the invariant mass distributions are fitted letting 
all parameters float, ii) the fits to the three angular distributions are performed fixing the Rj parameters to the 
results obtained in the previous fits. Hi) the invariant mass fits are repeated by fixing the 9p and ipp distributions 
to the results of the angular fits. Lpuj and 9^^ are always set to zero. 

Mass and widths of the p, w, and X are fixed. To account for the X width in the invariant mass distributions, 
we extract randomly the values of the mass of the X from a Breit-Wigner centered at mx = 3872 MeV with 
Tx = 1.7 MeV @, accepting only those values kinematically consistent with the decay of interest. 

Finally, since as described in Appendix |B] the invariant mass fits return consistent results with either the 
n = 1 or the n = 2 hypotheses (where n is defined in Eq. ([7])) we will only report results with n — 1. 

The resulting fits are shown in Figs. [5] and [3] and the fitted parameters are summarized in Tab. H] It is 
interesting to note that the radii, which are the fit parameters with a physical content, have reasonably small 
errors and get values consistent with 1 fm, the size scale of a standard hadron. Here it is hard to judge if we are 
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probing the size of the in which case we would say that the resuhs obtained for R disfavor a large 10 fm) 
loosely bound molecule [Q], or the size of the interaction making the X decay into J/ip + V. For sure the loosely 
bound molecule is requested to have a large wave function at the origin to make such a decay possible ^. On 
the other hand in any compact multiquark model the decay into J/tp does not require special conditions. 




0.55 0.60 0.65 

M(27i) (GeV) 



0.745 0.750 0.755 0.760 0.765 



0.770 0.775 0.780 0.785 



M(37l) (GeV) 



FIG. 2: Fit to the m^TT (left) and ms^ (right) distributions as described in the text, with the n 
curve refers to the 1"''^ hypothesis whereas the solid one is for the 2""*". 



1 model. The dashed 



A. Statistical interpretation 

The naive approach to evaluate the likelihood of the two spin hypotheses would be to consider the of the fits, 
obtained by adding the contributions from all the distributions, and to compare them with the number of degrees 
of freedom. In our case this would return a x^/DOF — 31.8/36 for the 1++ hypothesis and x^/DOF — 37.3/33 
for the 2 ^ hypothesis. The probability of obtaining a worse value is P(x^) = J°t P'DF{x, NuoF)dx = 67% 
and 28% respectively, i.e., both hypotheses seem to fit data well. 

Nonetheless this conclusion would not consider two effects. First, the overall good agreements is caused by 
the fact that we are simultaneously considering the distributions which favor the 1+^ hypothesis {9x, X: ^0 ^^'^ 
the one which favors the 2 ^ hypothesis (msTr)- This can be qualitatively seen in the fit results (Figs. [2] and E]) 
and in the Pix^) values obtained on each fitted distribution as listed in Tab.lTTl Then, a large number of degrees 





1^^ (dashed curve) 


2 ^ (solid curve) 


rjp 


0.089 ± 0.006 a.u. 


0.69 ±0.13 a.u. 




0.0026 ± 0.0003 a.u. 


0.030 ± 0.016 a.u. 


^ Jang 


1.32 ±0.4 a.u. 


1.03 ± 0.04 a.u. 






(254 ± 16)° 






(14 ±60)° 


R.J 


1.6 ±0.3 GeV-' 


5.6±0.8GeV-' 


xVdof 


31.8/36 


37.3/33 




67% 


28% 



TABLE I: Fit results for the two J assignments. 



Moreover we remind that loosely bound molecules are very difficult to be produced with high cross sections at hadron colliders [To| 
if not invoking final state interactions effect [TH . 
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FIG. 3: Fit to the angular distributions in X ^ J/tpTr'^Tv decays as described in the text. The dashed curve refers to 

the l^"*" hypothesis whereas the solid one is for the 2~^. 



of freedom are not sensitive to the spin of the X and therefore wash out the overall • This is the case for the 
m27r distribution, which is sensitive only to the radii, but not to the spin. 

To quantify this latter statement and develop a sounder statistical analysis, we performed a Monte Carlo 
(MC) study. We generate N data samples with the same number of events as the experimental samples and 
take into account the background, the statistical fluctuations, and the uncertainties in the model parameters. 
The simulations can be performed either by making the hypothesis that the X is truly a 1+"'" state or a 2 
state. 

As the starting point of the simulation for a given J and n hypothesis we take the corresponding model 
extracting sample by sample the parameters according to the result of the combined fit to the data. The 
invariant mass MC samples are generated by filling the bins bi according to Poisson distributions with mean 
values fii corresponding, bin per bin, to the sum of the values expected from the model and the experimental 
background. The angular pseudodata are generated as MC unweighted events. 

After the extraction this background is treated in the same way as the data, subtracted in the case of the 
invariant mass distributions and accounted for in the fits in the case of the angular distributions. 

The iV data samples obtained with this procedure are then analyzed with the same fitting function and 
statistical analysis as the experimental data. 

With such a tool, we first show the -P(x^) values obtained on the m27r fits (Fig.|3]) performed by generating the 
MC events with both the spin hypotheses and by fitting them either with the correct or the wrong hypothesis. 
These figures show that even when fitting with the incorrect X spin the fit would return a good P(x^). We 
therefore conclude that m27r is not sensitive to the X spin. 

Since we need the invariant mass distributions to perform the combined fit to constraint the Rj variables, we 
developed a more robust method of testing the hypotheses by using as estimator the difference between the 
obtained from the combined fits performed under the two J hypotheses [l^l 

Ax^=xMl++)-xM2-+) (26) 

Fig. [5] shows that the distribution of such a variable is peaked around a negative value when the MC samples 
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FIG. 4: Distribution of the P(x^) resulting from the fits to in MC samples generated assuming 1^"*" (left) and 2 
(right). The solid (dashed) histogram corresponds to the fit to the 2 ^ (l""^^) model. 



are generated with 1++ and, viceversa, it is peaked around a positive value when the samples are generated 
with 2 ^. These distributions are used to calculate the fraction of samples in which Ax^ has a value larger 
(for 1^"*") or smaller (for 2 ^) than the one obtained on data. This fraction, that we call CL, estimates the 
probability of the hypothesis. 

Given that in the combined fit Ax^ = —5.5, the 2 ^ hypothesis is excluded at 99.9% C.L., while the 1'''+ 
hypothesis has a C.L. of 5.5%. The different conclusion with respect to the naive expectations is due to the 
fact that individual values account only for the agreement of the data with the specific hypothesis under 
test. The Ax^ approach instead also considers the level of agreement with the other hypothesis, weighting the 
degrees of freedom with their power to discriminate among the two hypotheses. 



0.04 



0.03 

■a 

0.02 



0.01 



0.00 




FIG. 5: Distribution of the Ax^ = x^(l^^) — X^(2~^) resulting from the combined fits to Monte Carlo data samples 
with n — 1. The solid (dashed) histogram corresponds to events generated assuming the X to be a 2~^ (l"""^) state. 

We mark with a line the position of the experimental Ax^. 



B. Fits on sub-samples 

The Ax^ interpretation of the combined fit shows that even the agreement of the data with the 1++ hypothesis 
is marginal and the fit results seem to hint that the disagreement concentrates in the tostt distribution (see Fig. [5]). 
To quantify this statement we have performed the toy MC analysis on the Ax^ value obtained separately 
on the 7713^ distribution (obtained on the X ^ J/ip Tr+Tr^Tr" sample) and on the combination of the other 
distributions (obtained on the X ^ J/ip tt+tt^ sample). The results, listed in Tab. Illl show that the fits to the 
X — 7> J/ip 7r+7r~ distributions (second row) exclude the 2 ^ hypothesis and have an agreement with the 1++ 
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1 ~^ ^ ^/^QC?Vi^/^ /^ii"V"iro^ 

± lU-dbiifcJu. cuivt;) 




combined 


xVDOF = 31.8/36 
P (x') = 67% 
(7L — 5 5% 


xVDOF = 37.3/33 
P (x') = 28% 
(7 r — 1 % 


2tt (angular + mass) 


xVDOF = 20.9/31 
P (x') = 91% 

O-L = 23/0 


xVDOF = 34.7/29 
P (x') = 21% 

Oiv < U.i 70 


Stt (mass) 


xVDOF = 9.9/4 
P (x') = 4% 

O-Zj = U.i /o 


xVDOF = 1.5/3 
P (x') = 68% 
OL = oi/o 


combined (only mass) 


xVDOF = 25.2/22 
P (x') = 29% 
O — U . i /o 


xVDOF = 17.7/20 
P (x') = 61% 

U-L — 4D/o 


2n (only angular) 


xVDOF = 6.6/14 
P (x') = 95% 
Ci = 27% 


xVDOF = 19.6/12 
P (x') = 7.6% 
CL < 0.1% 



TABLE II: Results of the Toy MC 



hypotfiesis at 23% C.L., wfiile conversely tlie X — J> J/ip tt^tt^tt'^ data (third row) favor the 2 at 81% and 
exclude the 1+"'" hypothesis at 99.9% C.L. The two samples return therefore inconsistent answers. 

Finally, in order to compare directly with Ref. Q we report in Tab. also the results performed using 
exclusively the two invariant mass distributions or the three angular distributions. Also here the former give a 
clear indication in favor of the 2 ^ hypothesis, the latter in favor of the 1++ hypothesis. The different conclusion 
with respect to Ref. [Bl comes both from a different theoretical model and from the appropriate treatment of 
the degrees of freedom with no sensitivity to the spin of the X. 

V. CONCLUSIONS 

We re-analyzed the X — > J/i/j tt+tt^tt" and X ^ J /ij) TT~^Tr~ invariant mass and angle distributions published 
by the Belle and BaBar [Ij] collaborations respectively, with the goal to extract the most information about 
the spin of the X particle. 

With respect to the existing analyzes, we have improved two aspects. On one side, the X decay amplitudes are 
parameterized by effective strong couplings, weighting terms written as products of momenta and polarizations 
as dictated by Lorentz invariance, and parity considerations in a model independent way. The strong couplings 
are given a momentum dependency according to the model defined in Eq. ([7]). This requires the introduction 
of an additional parameter R which can be related to the finite size of hadrons in the interaction region. The 
results, shown in Tab. HI are consistent with R of the order of 1 fm. The treatment is fully relativistic and 
particularly appropriate for angular analyses. 

On the other side, in order to properly account for the sensitivities to the X spin of the considered distribu- 
tions, we pursue a statistical approach based on Monte Carlo simulations. 

We performed i) a global fit based on the whole information available from the 27r, Stt invariant mass spec- 
tra and the angular distributions of the X — >■ J/ip p decays, ii) two separated fits relative to the channels 
J/ip p (where we do the combined fit of the angular and invariant mass distributions) and J/ip uj (fitting only 
the invariant mass distribution). We also studied the fit to the invariant mass and the angular distributions 
separately. 

The combined fit i) excludes the 2 ^ hypothesis at 99.9% C.L., but returns a probability of only 5.5% of the 
1++ hypothesis being correct. The separate fits ii), return a clear a preference for the 1++ hypothesis in the 
J/ip p channel with a probability of 23% and an 81% preference for the 2 ^ assignment in the J/ip w channel. 
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FIG. 6: Feynman diagrams of X — >■ J/t/; p including the p-u} mixing. The second diagram includes the coupling e and 

the propagator of the u). 



with very strong exclusions of other hypotheses. Such results go in the direction of two different assignments 
for the two samples, which can lead to a host of interesting considerations. 
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Appendix A: p-to mixing 

As discussed in Sec. IIIII the isosping-violating p-ijj mixing introduces a correction near the pole of the w, 
without affecting the core of our analysis. Moreover, we found that the mixing improves the quality of the fit 
of 2 ^ and worsens the 

First of all, we have to insert the mixing into (|10p and ((T4)) . Following Ref. Q, we describe the decay uj — ?> 27r, 
as the oscillation uj p —?' 2tt, as explained in Fig. [S] We call — e the coupling of the vertex puj, whose value 
can be extracted by d, : 



m„ mp Fp F^ BTZ{uj 2tt) « 3.4 • 10"^ GeV^ (Al) 

The choice to treat the mixing reproduces naturally the phase of 95° of the complex mixing parameter used 
in Ref. 0]. 

Under the hypothesis 1++ we have 

{iIj p{q)\X) ^ gi^pT ~ egi^^— 2-7- (^2) 

where T is the S-wave scalar described in ([T]). Similarly for 2 ^ 

p{q)\X) = {g,,p + g^,p Ts) - e ^ ^1.^"'""^'' (^3) 

We can repeat the same argument for {ipu!{q)\X) by swapping the role of p and uj. The rest of Eqs. ([T3|) and 
(|17p (namely the Breit-Wigner and the decay into pions) is unchanged, because in this picture it is only the p 
(resp. the uj) which can decay in 27r (37r), being the mixing exhausted in the (ipVlX) matrix element. 

Moreover, if we take the mixing into account, the parameters rjp and rji^ of the fit appear together in both 
channels; therefore we must impose their correct relative normalization. This can be solved by imposing that 
T{X ^ J/tpSir) /r{X J/V'27r) = 0.8 ±0.3 

The angular distributions are unaffected by this term since they depend exclusively on the spin of the particles 
involved and not on the invariant mass spectrum. Therefore, to show the impact of the p-uj mixing we perform 
fits to the m2iT and distributions including this effect. The results are shown in Fig. [3 with the individual 
components detailed for the m2,r distribution in Fig. [51 The resulting values of the radii vary from Ri = 1.6±0.3 
to i?i = 1.1 ± 0.4 and from R2 = 5.6 ± 0.8 to R2 = 4.5 ± 0.7. As far as the hypothesis testing is concerned 
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the xVDOF changes from 25.2/22 to 29.3/22 for the 1++ hypothesis, and from 17.7/20 to 16.0/20 for the 2"+ 
one. None of the relevant conclusions is altered by this small effect and we will therefore neglect it in the main 
analysis. 




M(27i) (GeV) M(37i) (GeV) 



FIG. 7: Fit including p-uj mixing to the m2,r (left) and mz-n (right) distributions, with the model with n = 1. The 
dashed curve refers to the 1^"'" hypothesis, the solid one to the 2 ^ one. 

One could have thought that the presence of the narrow propagator of the lo would have constrained the fit 
to raise a peak at ttlui — 782 MeV. The smallness of e, instead, does not allow this, and the curve stays smooth 
at the pole, coherently with Refs. 0, H, 0. 




M(;t7i) (GeV) M(7t7i) (GeV) 

FIG. 8: Fit to the m27r distribution including p-ui mixing, with the model with n= 1, and under the hypothesis 
jPC ^ -^++ ^jgf^^ 2"+ (right). The solid light curve is the p contribution to the fit; the dotted light curve is the uj 
contribution, and the dashed curve is the interference term. The darker curve is the sum of all contributions. 



Appendix B: n = 2 case 

The theoretical model chosen for the fits described in this paper does not specify the value of n in Eq. ([7]), 
and it is therefore a free parameter of the model itself. Before adopting n = 1 for the rest of our considerations, 
we have performed a check on the dependence of the results on the choice of n. 

Since the angular distributions have shown a scarce dependence on _R, we fitted the invariant mass spectra 
with n = 2 comparing the result with the n — 1 case. The results are shown in Fig. [HI The resulting values of 
the radn vary from i?i = 1.6 ± 0.3 to Ri = 1.1 ± 0.2 and from i?2 = 5.6 ± 0.8 to R2 = 2.8 ± 0.3. As far as the 
hypothesis testing is concerned the x^/DOF changes from 25.2/22 to 24.8/22 for the 1++ hypothesis, and from 
17.7/20 to 18.9/20 for the 2 one. Also in this case none of the relevant conclusions is altered. 
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0.55 0.60 0.65 
M(27l) (GeV) 



> 10 

0) 



0.745 0.750 0.755 0.760 0.765 0.770 0.775 0.780 0.785 

M(37i) (GeV) 



FIG. 9: Fit to the m2,r (left) and 7713.^ (right) distributions as described in the text, with the model with n = 2. The 
dashed curve refers to the l"""^ hypothesis, the solid one is for the one. 
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